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SOUND RADIATION FROM LINE ARRAYS 


ieee Introduction 

The purpose of this paper is to develop simple approximations 
for the sound radiated from slender spheroidal bodies. The analysis 
of such bodies can be directly applied to the fundamental modes of 
submarine hulls and is therefore of practical interest in the field 
of naval architecture. Although the technique can be extended to 
any mode, this paper will deal only with the accordion and whipping 
modes, which are of the most practical interest. 

Although exact spheroidal radiation can be predicted from the 
available literature, the wave functions involved are unfamiliar and 
cumbersome. In addition to the inherent complication of the functions 
themselves, each author who deals with them tends to select his own 
notation, further exasperating the novice in their use. One signifi- 
cant result of this study has been my familiarity with these functions, 
however Slowly and painfully ence I have also gained an apprecia- 
tion of the desire to replace these functions by more convenient pre- 
dictions, which is the express purpose of this investigation. If one 
future investigator of such radiation is spared the use of spheroidal 
functions, this study will have fulfilled its purpose. However, since 
spheroidals form the exact answer that this paper is attempting to 
approximate, their use is mandatory to demonstrate the validity of 
those approximations. A brief summary of the necessary properties ou 


spheroidal wave functions is included in Appendix A. 


Literature on spheroidal functions is relatively complete. 

Gc. alerts and J. A. Senaceona have both published tables of 
spheroidal functions, but neither is ideal. Although containing 

all the required data, Flammer's tables are of rather limited range, 
and Stratton publishes only the expansion coefficients, requiring 
calculation of the other data from them. G. Bheeoe oe paper 
includes a concise and useful summary of the functions and analyzes 
some properties of individual spheroidal modes for rigid body axial 
vibration and accordion mode vibration. These analyses are reproduced 
ag ye sande eete” chaper on spheroidal functions, along with a complete 
and detailed development of spheroidal properties. A combination of 
all these sources was required to overcome the shortcomings of each 
individually and the variance in notation between them. 

As mentioned previously, this paper will develop approximations 
for the total sound radiation from the accordian and whipping modes of 
a slender spheroidal shell. Chertock's accordion results are for 
individual modes only, and not for the total radiation. M. eteacnerae 
has investigated the accordion mode radiation without Poisson effect, 
but demonstrates only fair agreement at higher frequencies. W. Sales. 
has analytically predicted the acoustic radiation from the transverse 


motion of a long, slender beam with elliptical cross section. Finally, 


C. Flammer, Spheroidal Wave Function (Stanford, 1957). 
J. A. Stratton, et al., Spheroidal Wave Functions (Cambridge, 1956). 
Ge 


i 

2 

3 ae 
Chertock, "Sound Radiation From Prolate Spheroids," JASA 33(1961) 871. 

4 

E. Skudrzyk, The Foundations of Acoustics (New York, 1971). 


M. Strasberg, "Sound Radiation from Slender Bodies in Axisymmetric 
Vibrations," Fourth International Congress on Acoustics, 1962. 


6 panes Jyh 
W. K. Blake, “Radiation from Free-Free Beams Under Influences of 


Light and Heavy Fluid Loading," NSRDC Report 3716 (1971). 





M. Teese” has formulated the basis for the whipping mode approximation 
developed in this study. This paper will extend the results of Strasberg 


and Junger and verify those results with the exact spheroidal radiation. 


7 eerie 
M. C. Junger, "Sound Radiation by Resonances of Free-Free Beams," 


MISA 52 (1972); 332, 





it. Procedure 

In general, the procedure followed by this paper will be to 
formulate an approximation of sound radiated from a spheroid by 
utilizing line arrays of sources with varying strength along the 
line. For the accordion mode, which is the longitudinal vibration of 
a body, the appropriate source type is the simple source; while 
for the whipping mode or tranverse vibration, a dipole source is 
applicable. Higher modes can be similarly approximated by line 
arrays of multipole sources. This approximation is then compared 
with its spheroidal equivalent. 

The first step of the approximation is definition of the shape 
and certain properties of the radiating body in both cartesian and 
spheroidal coordinates. Since the body used is a spheroid or an 
ellipsoid of revolution, this step is basic. The eee step is to 
formulate the incremental radiation from an element of length dz along 
the axis of the radiator. Such formulations are available from basic 
theory of sound sources. A surface velocity distribution is then 
selected that approximates the mode under investigation, and from 
it, the source strength is determined as a function of position on 
the axis. Finally, the incremental radiation is integrated along 
the length of the body, and an expression for the total radiation 
is obtained. 

The exact spheroidal radiation is determined as a summation of 
the individual spheroidal mode radiations. The acoustic pressure 
is related to the normal velocity, and the modal pressure coefficients 
are determined. Their evaluation is accomplished by matching the 


-velocity expression to the prescribed velocity distribution over the 


body and utilizing the orthogonality of the angle functions. The 
resulting integrals are then evaluated, and approximations for far-field 
and slender body are applied, resulting in expressions for each mode. 
These are summed to obtain total radiation. 

In principle, both formulations are straightforward but are 
complicated in practice by difficult integrations and infinite 
summations. Evaluation of the resulting expressions is also cumbersome, 
involving extensive numerical manipulations. As will be pointed out 
later, some effort must be devoted to maintain precision, if ficticious 
peaks in the smooth functions are to be avoided. 

Since the two methods are obviously considerably different, the 
arrival at equivalent expressions fortunately is an indication that 


both are probably correct. 





III. Accordion Mode Approximation 

As pointed out in the procedure section, the first step in the 
formulation of the problem is to develop certain properties of the 
radiating body, in this case, a spheroid. The expressions developed 
will be used in both the accordion and whipping mode approximations. 

For an ellipse of focal length dad, major axis L, and minor axis 
D, a2 = L4 (1 - (D/L) +), as shown in Appendix A. Therefore, for 
slender spheroids (D << L), dis approximately equal to L. For 
purposes of this discussion, the other properties of interest can be 
derived for the ellipse that is the intersection of the spheroid and 
the x-z plane. In this case, x is the perpendicular distance to the 
axis and z is the distance along the axis with z = O at the center 


of the body. The governing equation for this ellipse is: 


Z Me 
2 eee 2 al 


(L/2)* (D/2) 
If a(z) is defined as the radius of the cross-sectional area 


normal to the axis, then: 


We 


a(z) = x(z) = DIZ - (2/L)) (1) 


If a(z) is the angle between the normal to the spheroid and 


the normal to the z axis, then: 


tan(az eaten /d2e = (2) 


L2 (= 2 (Gym 29 
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Using that definition of a(z), the normal velocity of the 
spheroidal surface is made up of contributions from both the axial 


and "radial" (perpendicular to axial) velocities as follows; 


ue 7 u,sin(a(z)) + u cos (a (2) ) (3) 


For the array approximation, only the tangent form is required; 
however, for the spheroidal work, the individual sine and cosine 


terms must be obtained. Using: 


sin(a(z)) = me 
USE GQ 
cos(a(z)) = ox 


(lax? + (ae)2)—7— 


and differentials obtained from Skudrzyk's aotaeaione® leads to 


the following relationships for the spheroidal velocities: 


sin(a(n)) (4) 


i 
= 





Ay hye O 


cos (a(n) ) (L/D) (1 -n 


(5) 


Finally, the incremental surface area along length dz is given 


by: 
dA = 2Ta(z) ds = 2tTa(z) dz/cos(a(z) ) (6) 
The radiation from a simple source is simply: 


p = (p exp(ikr)/4mr) § (7) 
jr=a 


2s ugha, perdog,,eq, 14, 15, and 16. 
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where p is the acoustic pressure, p is the ambient density, k is the 

wave number, r is the distance from the source to the point of 
observation, and S is the volume velocity (S = dS/dt). The incremental 
volume velocity for an element of length dz is given by ds (z)=u_ (z) daA(z), 
and using (3), (6), and (7), the incremental pressure for that element 

is: 


dp = (-ipwexp (ikr)/2r) a(z) [ue tu, (2) tan (a(2))] ae (8) 


adopting exp(-iwt) time dependence. 

For the far field approximation, small differences in range 
are ignored in the magnitude term, but retained in the phase. Hence, 
if r = R- zcos8, where R (R >> L) is the distance from the center 
of the body and 6 is the angle with the z axis, the total pressure 


is expressed as: 


+L/2 

p =<.(1pwexp (1kR)D/2R) f exp- (ikz cos (8)) (> - (eye 
-L/2 
u_ (2) Dez 

X Cu (z) 7 ee) daz (9) 


L2 (5 - (2/L)2)7/* 


where use has been made of (1), (2), and (8). 

The accordiOn mode of vibration is the fundamental, axial, 
elastic mode, with maximum velocity at the ends, no axial velocity 
at the center, and a lateral ("radial") velocity due to Poisson 
effect. This lateral velocity is maximum at the center. For a 
uniform tube, the ratio of lateral to axial velocity is -ToD/2L, 
where 0 is Poisson's ratio. Therefore, for a thin spheroidal shell, a 
fair approximation to the motion of the accordion mode is: 


=| Shercasrds, Dao ron 


a2 





u, (2) Uy Sin (1z/L) (10a) 


uy (2) ~(ToD/2L) Up cos (12/L) (10b) 


Substituting these expressions into (9) yields the following 


expression for non-dimensional pressure P , 
fe) 


P ee (11) 
(exp (1kR) (-iwU,) (D/2)? 
where: 
+L/2 
fi = one f exp(—1kz2 cos()) 2 sins (7z/b) az 
-L/2 
-L/2 
- iL eae 
In = -10/L (7 > (Z7 i} exp(-1kz cos(8)) cos(mz/L) dz 
~-L/2 


: ; JAE 
Both integrals can be evaluated in closed form, thereby ele mee 


I, = (w2-v2) 72 [wr tw2—v2) sin (w)+(w2+v2) cos (w)] Ci) 


J) (w-v) J) (w+v) 


In = (-14%0/8) Sama ie os ~ Garey) 


where w = kL cos (8)/2, v = 1/2, and J; is the Bessel function of the 
first kind, of order 1. I), agrees exactly with Strasberg's result, 
but he has assumed o = 0, so that u = O, Ip = O for his analysis. 
Combining (10), (11), and (12) leads to the desired expression for 


the total sound radiated by a spheroid in accordion mode vibrations; 


a. SCE, p. 968). eq. o>. 


css coerasansi fer peo oOmmca. 9.1. 20., and p. 225, eq.6o6.1.9'. 
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Py = (w2-v2) ~2 { w (wv?) sin(w)+(w2+v2) cos (w)] 


J} (w-v) J) (w+v) 


m2 On eee 
(r%o/8) ( (w-v) iu (w+v) 


Evaluation of this expression gives the smooth curve on Figures 
1 and 2. In numerical calculation, care must be exercised to maintain 
precision in the region where w approaches v (w = v + €), Since terms 


of order e2 


must be retained within the parentheses to avoid a ficticious 
peak in the curve. If precision is maintained, the first expression 
converges to a finite value for w = v. The curves for I] and Ip are 
plotted separately in the graph to allow individual comparison with 
spheroidal predictions and Strasberg's results. 

For calculation of I5, a Poisson's ratio of .29 was used, which 
is the static value for a uniform steel shell. However, in practice, 
the value used would probably be modified by the addition of stiffeners 
to the shell and by dynamic effects of ring resonances at higher 


frequencies. Since the same value was used in both calculations, 


these modifications have no effect on the results here. 
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IV. Accordion Mode Spheroidal Calculation 

The spheroidal wave function equivalent of the preceding expression 
will now be evaluated using the spheroidal functions that are outlined 
in Appendix A. For purposes of brevity, the abbreviated symbols 


a and Ron will be used to indicate the complete variables San en) 


and R BP rn, , Rnle will be used to indicate ae ie) 
0 


on the surface of the radiator. 


(h,&) evaluated 


The total acoustic pressure is the sum of the pressure contri- 
butions from each of the individual modes. If Pian are the modal 


spheroidal pressure coefficients, then: 
“) ‘y Pan Xmn °mn °°S (mo) (13) 
m n 


The relationship between pressure and the normal component of 


particle velocity is given eves 


where Fs is the metric sea BRi Asie, Mnesetone: 
Ld — Renemn COS im 9) = iweu, (n, >) Te (14) 


where Ron = OR! o8: If this equation is applied as a boundary 
condition over the surface of the radiator, the orthogonality of 
the angle functions can be utilized. Multiplying both sides of (14), 


nes onsbaepay p. 467. 


Pol sienabaaele, D. 459,060. LS. 


ILS) 





evaluated on BE oe by Su cos(k¢) and integrating over the surface 


yields: 


27n +1 


D2. aie J J cos (md) S,, cos(ko) dn dd 
re) 


27 +1 


= iwo { ] a pwns S44 cos (ko) dn dd (ES) 
E 


‘e) 


The dnlepeal on the left hand side of the above equation reduces 
to a multiple of the normalization constant and vanishes unless 
=k andn= 1. The only non-zero value is 27™N /e where € = 1 
Tid, mi’ Oo 
and ii 2 for m # Q, and Nan is given by equation (A.3). 


Combining (13) and (15) and making use of (A.2) leads to the 


following expression: 








cos (md) 
p = caus wae SB i 
Aa Ron 
Bs 
2n +1 
ae 72 =e 
am f J UL Peet) cos(m¢)dndd (16) 
s 1 55 | 
LS 


which agrees with Chertock. 
For the accordiOn mode, the spheroidal modes that contribute to 
the radiation are the m = 0 and n = even modes, since such vibration 


“dchevabeayey p. 470, eq. 56. 


5 
: Ghemtock,, pova’s, eq. zc and 13. 


16 





is independent of ¢ and symmetric about the midsection. Combining 
(16), (A.4), (A.5a), and (A.5b) with long, thin spheroid (d>*L) and 


far field (&+2R/L,n~>cos8) approximations and integrating over ¢ leads to: 


co 


2 : On,.~-(n+1) n TS 
p = (woL*exp (ikR) /4R) > (do (2. is n° ((n/2) 1) (Es 1)} /No nt 





O 
n=0 
even 
eo +1 
) on if 2 ai We 
X d ———-—- Wa (hp Pa) an (17) 
Xr a f 2a ai r 
r=0 fe) fe) 
even 


An expression for the normal velocity on the surface can be 
derived from (3), (4), and (5) in spheroidal coordinates and is 


given by: 


Ae 
Bo 71 


yo 
u_(n) = U. |- ( sin(m/2) -(10/2) (1-n?) 
n O E ae 

O 


eos tan/2 


(18) 
which agrees with enertoceecs Combining (17) and (18) and using the 


earlier definition of P (equation (11)), then: 


=a One os. n 2 
Po = (do sin((n+)1/2)S, 2 COnr 2 el) /N oy”! 
=Q 
even 
On 
4 Sag ae) (19) 
r=0 
even 


MC arahocts: Wee ovOy CG. 632: 


17 





where: 


+1 
nas if nsin(mn/2) P(n) dn 
-1 
+1 
as | (ee n2)*/4 cos(m/2) P(n) dn 


~l 


The integral terms J. and KE were then evaluated for r = 0,2,4 
by expansion of the Legendre polynomials and doing the resultant 
groups of integrals. Both groups are integrable in closed form, but 
the K series results in Bessel functions of the first kind of order 
r/2 + 1, which makes their evaluation somewhat cumbersome. The 


resultant values are: 


1g J K 

1g 1g 
0 2ol057 3 37. 
Z 28567 -,2970 
4 - .02763 ,0376 


Utilizing the above values and tabulated data on the spheroidal 
functions, equation (19) can be used to determine Poe It should 
be noted that the expansion coefficients (ao") in (19) are those 
given by Flammer. If the wider range of Stratton's tables is desired, 
the difference in normalization’ ’ must be accounted for. For use with 
Stratton's tables, equation (19) becomes: 


eoanness De 6G. 
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oO 4 


dy (alon) ) d_(n[On) P_(n) y d(h[On) (J - (10/2) K,) 


a = r=0 r=0 so 
n=0 oon men 
= zZ 
Gs yy 2(d_ (hon) )?/(2xr +1) 
r=0 
= = Wee x y) 
ao = x 1) Yr] d_(h|On) /2 (Crt) 
r=0 


In this case, Stratton's tables were used with the above equation 
and excellent agreement with the approximation was achieved (as shown 
by the crosses on figures 1 and 2 and as demonstrated by Table I.) 
Calculations using Flammer's tables and equation (19) also demonstrated 


the approximation‘’s accuracy. 
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FIGURE 2 
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TABLE I Accordion Mode Results 
i. Tha Ths 
0 ~4052 - 4053 
220 .4007 -4007 
5218) RSIS: - 3769 
A 1sh0) 3343 . 3343 
1 .00 ~ 2966 . 2928 
JL alate, wld seb swales) 
fae PAW, -.0143 -.0169 
yt ehe, -.1619 -.1620 
3.40 -.2596 -,2591 
4 .40 -.2896 -.2895 
ir - First term of array approximation 
L1s - First term of spheroidal prediction 
I,4! I4, - Second terms of array and spheroidal, respectively. 
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V, Whipping Mode Approximation 
The whipping mode is the beamlike transverse vibration of 
the shell. Since a slender spheroidal shape is again assumed, 
the properties derived in section III are applicable. The radiation 
from an element of length dz of the spheroid can be developed 
from the expression for a circular cylinder given by Junger and 


on alts! 
Feit. For such an element, the n = 1 term is: 


ue) 
I 


(pa exp (ikR) /R) (2W, E(k cos (6)) (-1) (ka siné) cosd/4 


where: 


Z4+dZ 


£(k cosé@) = if cos(kz cos8@) dz = dz 
Z 


Therefore: 


dp = -wu_(z) pm (a(z)) “exp (ikR) ‘sin@ cosh dz/RA (20) 


which agrees with eneeen 


A transverse velocity distribution of the form U) (cos(1z/L)-2/t1) 
: , , 20 
would be the simplest vibratory mode of a uniform free-free beam, 
For a slender spheroidal thin shell, very little error is introduced 


by using that velocity distribution. Substituting into (20) and 


using (1) and the phase difference term leads to: 


ener ange helt, pe Zio. 


eerie aa Dae os2. 


20 nen Harte, jp. loa. 
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Ty 


Is 


I¢ 


5 Ree (ig + Ty + Is + ee. 


wU Lp 1Dexp (ikR) cos ($) 
L/2 
(1/4L) J cos(mz/L) exp(-1kz cosd¢) dz 
-L/2 
L/2 


SVEe) [ ereostrer exp(-1ikz cos@) dz 
-L/2 


L/2 
-(1/2TL) J exp (-ikz cos@) dz 
-L/2 


by 2 


(2/mL3) J z*exp (-ikz cos§@) dz 
Bye 


All four integrals can be Avra eens in closed form, yielding: 


13 


Ty 


I¢ 


= mtcos(w)/8(v2 - w2) 
= 81 [eos (w) (6w2 + 2v4 - (v2 - w2) 2) 


-~4w sin(w) (v4 - a (v* - w*)73 


- sin(w) /21w 


I! 


(w cos(w) + sin(w) (w2/2 - 1)) /mw? 


2 
peyadeneeray p. 198, eq. 8. 
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Combining the above expressions leads to the desired approximation: 


. 


Py = [way (v2 ~ w2)73 (cos (w) (3w2 + v2) = 2w sin(w) (v* - w7)) 


- sin(w)/2tw + (w cos(w) + sin(w) (w*/2 - 1)) /nw9 Sing 


Numerical evaluation of this expression for two different 
frequencies is presented as Figures 3 and 4. Once again, care 
must be used to maintain precision in the region where w approaches 
v. In this case, the numerator of the first term must be evaluated 


to terms of order e* to avoid a ficticious Tea ey. 
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VI. Whipping Mode Spheroidal Calculation 

Calculation of the radiation from whipping mode vibration using 
spheroidal wave functions makes use of the basic relationships developed 
in section IV. Equation (16) is the general expression for all modes 
which is specialized to the whipping mode. 

For the assumed whipping velocity distribution, the spheroidal 
modes contributing to the radiation are the m = 1 and symmetric 
modes. For m= 1, the symmetric modes are for n - m even, or ll, 13, 
etc. The complete expression for the normal velocity on the surface 


is given by: 


u_(2,9) = u (2) cos(a(z)) cos®¢ 


By making use of (5) and shifting to spheroidal coordinates, 


that expression becomes: 


aS DZ 
bo 72 


172 


u(n.¢) = (LU,/D) (cos (mn/2) - 2/n) (1-n*) ’ “cos¢ (21) 


Ded 
50 n 


Gomoaning (21) "and (16) with “(A.4), (A.5a), (A.5c), (A.5d) and 
long, thin spheroid (d*L) and far field (&>2R/L, n-cos8@) and integrating 


over ¢ leads to: 


~ go oaell), n yn 
4 -i S12 ((n-1)/2) 1 ((nt1) /2) 'd22B,_ A} 
P, = -(1/4h?) ) Premed Ns a ae 


Ny (n + 1)'Cy, 
n=1 


where: 


26 





co | andl 


Be -) aw? J xn22/? (cos teny2)-a/m?, dy (n) dn 


r=0 -1 


co veil 


Z in 2) - d 
al if (1-n2) (cos (1/2) -2/7) Gn (Prey (M)) dn 


r=0 -l 
co ! 
_ 1m 
B) = (r+2) (r+1l)d 
n Yr 
r=0 
0 I 
C). = ae 
n r 
r=0 


which has been non-dimensionalized by multiplication by l = ThL/ dh 


and use has been made of the Pediaetenensint << 


1 Oly 2a 
Pp — — ie 
r+] 2 Nee dan Prtl 
The integrals in Als can be evaluated using the same method used 
in the accordian mode section: expanding the Legendre polynomials 


and integrating. The values from those operations are: 


TsO, e322 
r= 2, -.490870 


r= 4, a hOPS MNS) 2h, 


In this case, the pressure was evaluated using Flammer's tables 
over the same range of frequency. Figures 3 and 4 are marked with crosses 


indicating the spheroidal prediction from h = 2 and 4 over a quadrant 


Deer: 
2 Haldebrand,p.. 165, eq. L72. 


A 





of @ variation. As with the accordiOn mode, the line array correlates 
very closely with the spheroidal results. Table II gives numerical 


results for comparison. 
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TABLE II - Whipping Mode Results 


0 P) (h=2) P) (h=2) P} (h=4) P) (h=4) 
Approx Exact Approx Exact 

15 .00720 00726 s007 24 .00739 
30 paotso 2 .OU356 -01466 .01484 
45 -01836 201836 gO2112 sO2) 
60 mo ney wOztZ7 2 430 .02417 
75 g02256 502257 ,02383 .02384 
90 702290 nO'2 290 202290 ~G2306 
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VII Near Field Approximations 

The most severe restriction assumed in the derivation of these 
expressions is the far-field approximations. These approximations 
simplify the evaluation of the expressions and permit solution in 
closed form. However, practical considerations encourage the develop- 
ment of models that can be used at observation locations closer to 
the radiator. Therefore, the evaluation of the approximations was 
extended to the near field. Time constraints permitted such evalua- 
tion for beam aspectath = 2 only. However, the computer program 
developed has purposely been kept general and can be used for any 
position or frequency. 

In the prior derivations, the far-field assumption permitted 
extraction outside the integral of the exponential and range terms 
(Equation (8) to equation (9) and application of equation (20)). 
Retaining these terms within the integral leads to integrals that are 
practically impossible to evaluate in closed form. Therefore, 
numerical evaluation was selected, and a computer program was developed. 


The integrals so evaluated are; 


+1/2 
Ee = DIST { axe iz sin(1Z) - — (1/4 - Boye co (1 az 
-1/2 


for the accordion mode, where Z = Z/L and DIST = R/L; and 


+1/2 
P, = DIST I Sos (2 (= - 2°) dz 
2 
-1/2 5 


a2 





for the whipping mode, where RF = Distance Abeam/Length. In both 
cases R must be retained as a function of Z for near field results. 
The computer program calculates and sums the real and imaginary 
contributions from each segment of length aZ along the body and 
then calculates the magnitude of the pressure term. A listing of 
the program is included in Appendix B. In all cases, convergence 
to less than 1.9% was required, but in most calculations, the con- 
vergence was considerably better. For all whipping calculations, 
convergence was less than 0.04%. 

The spheroidal values for the near field pressure magnitudes 
were calculated using the following expressions, which were derived 


from equations (19) and (22) for beam aspect (@ = 90°); 


By h DIST [227949 ~- .168R\5))? = (.127R'2) - sen?) )2| we 


00 2 00 02 


ty 
— 
| 


a (1) CR) (2) (2), 51 eZ 
DIST/h [(-206R6} + .022R,,°)* + (.206R,)° + .022R,, 4 


Values for the radial spheroidal functions were obtained from the 
NRL es which give only limited data for values of &greater 
than 1.4. For both modes, only the first two spheroidal harmonics 
were used: 00 and 02 for accordian and 11 and 13 for whipping. 
The values calculated for both the approximations (solid 
curve) and spheroidal expressions (crosses) are shown in Figures 
> and 6 and Table III. The far field computer results for both modes 
agree with solutions of prior sections. The agreement for the accordion 


mode is extremely good to a distance of a quarter-length abeam. The 


23 : 
Hanish, S., et al., Tables of Radial Spheroidal Wave Functions, 


NRL Report 7088 (Washington, 1970). 
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whipping mode is not as good and can be useful to a minimum distance 
of about half-length. Since convergence for whipping numerical 
calculations is very good, this error can be attributed to a 
combination of inherent inaccuracy of the model and the possible 
inaccuracy introduced by truncating the spheroidal evaluation at 
only two terms. Therefore, in view of the possible inaccuracies of 
each solution individually and the variation between them, neither 
near field whipping mode solution should be considered the exact or 


correct solution. 
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TABLE III - Near Field Results 


Distance Spheroidal Approximation 
Length (Convergence (%) ) 


Accordion Mode 


ALP LOFG, . 1446 L453". 65) 
Ag ie, i326 sooo. OO) 
529) wt 230 peso (240) 


Whipping Mode 


2.00 .02344 702325701) 
med, BO2519 502422(,01) 
75 502579 .02452(.001) 
Bo2 NOPKSS TAL 202489001) 
moO sO2e622 702523 02) 


oy, 


BYrror (s) 


48 


oe 


273 


ashy, 


5.00 


10.40 





VIII. Summary and Conclusions 

This paper has developed an approximation for the sound 
radiation of two fundamental vibratory modes of a spheroidal 
body. This approximation has demonstrated very good correlation 
with the exact spheroidal prediction over the frequency range 
of practical interest. The expressions developed are far more 
Straightforward than the cumbersome spheroidal wave functions. 
Their future use should therefore promote faster and easier 
predictions of such radiation. 

The two obvious limitations of these approximations in 
practical application are the near field performance and the 
physical limitation imposed by a spheroidal shape. Section VII 
has demonstrated excellent near field performance for the 
accordion mode approximation and fairly good performance for the 
whipping mode, at least in the intermediate or "suburban" field. 
The second limitation of shape differences between a spheroid 
and actual physical radiators should not be severe, particularly 
Since the actual shape can be incorporated into the approximation 


and evaluation carried out numerically. 
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Appendix A - Spheroidal Wave Functions 

Spheroidal wave functions are the result of the solution of 
the wave equation in spheroidal coordinates by the usual technique 
of separation of variables. The prolate spheroidal coordinates that 
are of interest here are obtained by rotating elliptic coordinates 
about the major axis. If the foci of the ellipse are located at 
iGaey, 2) = (0, 0, td/2) and rj and ro are the distances from these 


foci, then the spheroidal coordinates are: 


5 Sse /-e! 1<é<@ 
N= (rj-ro)/d STs 
¢ = tan7!(y/x) O<$<2T 


Surfaces of constant € are confocal ellipsoids; surfaces of 
constant n are confocal hyperboloids orthogonal to the ellipsoids. 


For a spheroid of length L and minor diameter D; 
e* = 1 + D*/a? 


ane liipse issdetiined by Surfaces where rn) toro = Jconstant. 
Hemece, On the major axis, rr; + ro = L, and on the minor axis, Fr), + F929 = 


2((D/2)2 + (a/2)2) 142, 


Therefore, L? = Dp? + 2, andl FOR, ->. Deer 
for a slender or needle shaped spheroid, d> L. The radial coordinate 
on the surface, or EO! approaches 1 for slender spheroids. 

Pn ENee tateee ld) esq 690 sea ee ey lpm whe ee tvert Seen 


average distance to the spheroid. Furthermore, in the far field, 


Nn > cos (8), where @ is the angle with the major axis. 
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The separable solution is of the form: 


aa S an en) Ronee) cos (md) (A.1) 


where the indices m and n (m <n) identify the various modes. The 
parameter h is a non-dimensional frequency, h = Swd/c = Skd, where 

W and c are the angular frequency and sound speed, and k is the wave 
number, or their ratio. A surface harmonic is represented by 

. San en) cos (md) which vanishes at n-m circles of latitude. 

Through application of the Laplacian operator in spheroidal 


coordinates to the above solution, it can be shown that both a and 


en satisfy the same ordinary differential equation. 





d Qa SLL 2.2 ue fo 
[v= S| j Inez Jere 


mere tor -l1<z<+tTI, Z=yj and u= §S 


ema cor tl < Z < &, z= € and u = Ron 


The angle functions San (hen) are generalizations of Legendre 
functions and depend onm,n,h, and n, thereby making spheroidal 
Manipulations cumbersome. They are generally expressed as series 


expansions of Legendre functions as follows: 


ao if 


mn m 
San fen) = 3 ele) ee are) (A. 2) 
x=0,1 


where the prime over the summation sign indicates summation over only 


even values of r when n-m is even, and over only odd values of r when 
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n-m is odd. The expansion coefficients are available in Flammer and 
Stratton, although the range of the latter is much greater. 
A quantity required for calculations is the norm or mean square, 


which is also a function of m,n, and h, and is given by: 


+1 co | 
f : (r+2m) | [an (h) }2 

Nan (?? = (Shen) dn = 2 Moo ae (A.3) 
-1 r=0,1 


This quantity is tabulated for limited ranges in Flammer, but must 
be calculated from the expansion coefficients in the use of Stratton's 
tables. 

The radial function is a complex function describing the distance 
dependence of the wave functions. It may be computed at any point by 
a series expansion of spherical Hankel functions of the first kind. 
Additionally, numerous tables are in existence for a wide range of 
frequencies and distances. For time dependence exp(-iwt), which is 


used here, the radial solutions are: 


(3) Dees es) ue) 
a (h,€) = ISS iy et ae (hiya? 


which correspond to the diverging wave. 
However, for purposes of this study, only the far field approxima- 
tions for long and thin spheroids are required. These are available 


nm ooth hegate and Seer eee Ips, ee Le 


2p. Skudrzyk, p. 464, p. 475-476. 


41 





Ri.) (h,6) > y7 692) xp (GhE) / (hE) (A. 4) 


Since h = = kd and — > 2R/d in the far field, the exponential term 


becomes exp(ikR). 


A further quantity required for the calculation of sound radiation 


is the derivative of the radial function. For slender spheroids: 


' zg Pao). 
Ri.) (hy) = Re eo) > ir’) (n,é) BIS) 
eae 
For m = Q; 
(2 1(€2-1)°} 
Ron (hi &) > sa AULA ee ee n even (A.5(b)) 


meee? (ny72) (a /2yen a>" (nh) 


For m > Q: 


(2) = =1/2m- 

RA? (hE) + =m£g(myn,h) (€2-1)77/ 2 (a.5(c)) 

Aol 
Where: 
m n-1 mn 
(-1) (2m—1) (m-1) ! (m+1) th 3 d_. (h) 

fo(m,n,h) = ee n even 

jn-2mtl (2m) | cS Can 1mm we (h) y fomee) ! a” (hn) 


(A.5(d) ) 
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Appendix C 


a(z) 


eee Ul 
m 2 LY 


- Important Symbols 


radius of cross section 

sound speed in medium 

focal length 

non-dimensional frequency = ska 
wave number = 27/A 

index 

acoustic pressure 

time 

velocity components 

1/2 

kL cos (6) /2 

maximum diameter of spheroid 
length of spheroid 
normalization constant 

Legendre functions 

Accordion mode non-dimensional pressure 
Whipping mode non-dimensional pressure 
Far field range 

Radial spheroidal function 
Angle spheroidal function 
Maximum velocity, accordion 
Maximum velocity, whipping 
Radial coordinate 

Angle coordinate 

Longitudinal angle 


Angle with axis 


46 





Appendix D - Bibliography 


The followingsreterences were Consulted: 


ie 


iLO 


V1 


2. 


ILS 


14. 


TO. 


Abramowitz, M. and Stegun,I.A., ed., Handbook of Mathematical 
Functions (National Bureau of Standards, 1970). 


Blake, W. K., "Radiation From Free-Free Beams Under Influences 
of Light and Heavy Fluid Loading," NSRDC Report 3716(1971). 


Chertock, G., "Sound Radiation from Prolate Spheroids," JASA 33(1961). 
Den Hartog, J. P., Mechanical Vibrations (New York, 1956). 
Flammer, C., Spheroidal Wave Functions (Stanford, 1957). 


Gradshteyn, I. S. and Ryzhik, I. M., Table of Integrals, Series 
and Products (New York, 1965). 


Hanish, S., et al., Tables of Radial Spheroidal Wave Functions, 
NRL Report 7088 (Washington, 1970). 


Hildebrand, F. B., Advanced Calculus for Applications (Englewood 
CHicecetimmew, LIOZ). 


Junger, M. C., “Sound Radiation by Resonances of Free-Free Beams," 
JASA 52° 4@ig7 2) . 


Junger, M. C. and Feit, D., Sound, Structures, and Their Interaction 
(Cambridge, Mass., 1972). 


Morse, P. M. and Feshbach, H., Methods of Theoretical Physics 
(New York, 1953). 


Morse, P. M. and Ingard, K. U., Theoretical Acoustics (New York, 
1968). 


Skudrzyk, E., The Foundations of Acoustics (New York, 1971). 


Strasberg, M., “Sound Radiation from Slender Bodies in Axisymmetric 
Vibrations," Fourth International Congress on Acoustics, 1962. 


stratton, J. A., et al., Spheroidal Wave Functions (Cambridge, 
Mass., 1956). 


47 











Thesis 145674 
A3345 Akst 
Sound radiation from 
line arrays. 


thesA3345 
Sound radiation from line arrays. 


Ml il | i iit 





